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Abstract 

The conditions for zero-energy Andreev surface bound states to 
exist are found for the lattice model of d-wave superconductor with 
arbitrary surface orientation. Both nearest neighbors and next nearest 
neighbors models are considered. It is shown that the results are 
very sensitive to the surface orientation. In particular, for half-filled 
(/i/0)-surface zero-energy Andreev surface states only appear under 
the condition that h and I are odd simultaneously. 

Significant feature of high temperature superconductors (HTS) are zero- 
energy Andreev surface bound states. The zero-energy states (ZES) form 
on surfaces of a <i-wave superconductor with orientations different from the 
(100), due to the sign change of the order parameter. In high-temperature 
superconductors, such states manifest themselves as the zero-bias conduc- 
tance peak in tunneling spectroscopy in the ab-plane 01213111001111313 

cnj en nil iini eh cni cni im im iei eoi eh E21 ehi eh ehi eei , the anoma- 
lous temperature behavior of the Josephson critical current [23 [23 123 ED] 
and the upturn in the temperature dependence of the magnetic penetration 
depth [^TJ E2 ES] • At the same time the problem of ZES at the surface with 
arbitrary orientation is still not clear. 
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The conventional description of Andreev surface bound states, as well as 
the most of inhomogeneous superconducting problems, is based usually on 
the continuous quasiclassical approximation. From the this viewpoint the 
conditions for ZES to exist are quite simple. ZES are formed due changing of 
order parameter sign along the quasiclassical trajectory. There are no ZES 
for (100) (i.e. 0°) orientation, and there are ZES at all values ku for (110) 
(45°) orientation. For intermediate surface orientations the sign change does 
not take place for all incoming momentum directions, and the weight of the 
ZES decreases with deviation from 45°-orientation. 

From the other side, the tight-binding BCS model is widely used for 
theoretical description of HTS. This model gives the same (as the continuous 
quasiclassical model) result for ZES at (100) and (110) orientations. However, 
for the intermediate surface orientations the question is very complicated. To 
the best of my knowledge only the simplest orientations (100), (110), (210) 
[2*H |23 OH HHH EH IHZj have been studied. And even numerical calculation 
couldn't give the general answer for all surface orientations because of the 
lattice specific. 

In this paper the general analytical criterion for zero-energy Andreev 
surface bound states to exist at the surface of arbitrary orientation is pre- 
sented. I consider two-dimensional tight-binding model on square lattice. 
The surface orientation is assumed to be arbitrary and characterized by the 
indexes (W0). Both nearest neighbors and next nearest neighbors models 
are considered. For simplicity I take the superconducting order parameter 
to be spatially constant. This approximation is reasonable for studying of 
low-energy quasiparticle states. The impenetrable surface is assumed to be 
smooth. 

The Hamiltonian for a pure singlet superconductor in the tight-binding 
model can be written as: 

n = - *(x,x')4(x) Cff (x') + 

x,x',cr 

+ ^{A(x,x04(x) C {(x')+/ i .c.}, (1) 

x,x' 

Here i(x, x) = \x is the chemical potential; £(x, x ± a) = t(x, x ± b) = t > 0, 
£(x, x±a±b) = t'<0 are the hopping elements; ci-wave superconducting 
pairing is defined for nearest neighbors A(x, x ± a) = — A(x, x ± b) = A. 
Here x - are the positions of lattice sites; a, b - are the basic lattice vectors. 
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Then the Bogoliubov-de Gennes equations take the form: 



-t(x,x') A(x,x') 
A*(x,x') t(x,x') 



We define new coordinates (x, y), rotated with respect to the crystal 
axes (a, b), where x is the direction normal to the surface and y is the 
direction along the surface. Superconductor is situated at x > 0. Lattice 
constant is taken to be unity, a = 1. The system is periodic along the y- 
direction with period \Jh 2 + I 2 = d^ 1 and the crystal momentum component 
k y of a quasiparticle is conserved. Instead of the usual square Brillouin zone 
(BZ) k a G [— 7t, 7r], k b G [— n, it] we now use the surface-adapted Brillouin 
zone (SABZ) (221 EH] given by k x G [—ir/d,ir/d] and k y G [—7id,nd]. Here 
d — 1/ \Jh? + I 2 is the distance between the nearest chains (layers) aligned 
along the surfaces, i.e. all x-coordinates have discrete values with period 
d. The momenta in the two coordinate systems are simply related through 
rotation of an ang le 9 = tan" 1 h /I. 

Let us solve Eq.([2|) for half-space x > and fixed k y . General solution is 
constructed from all the solutions of uniform problem, which don't grow at 
x — > +oo. Then the wave function for fixed k y can be written as 



K V(x,ky)J ^ \V a (ky) 

here summation should be taken over all solutions k x a of equation 

E 2 =e(k X ,ky) +A 2 (k X ,ky), (4) 

with Imk xa > 0. The boundary conditions are 



u(-jd, k y ] 
v(-jd, k y ] 



j=0,l,...,iV-l, , (5) 



where = max(/i, I) for nearest neighbors model (t ^ 0, t' = 0) or N = 
h + I for next nearest neighbors model (t, f ^ 0) [22] • The total number 
of solutions (jlj) with lmk x a > equals to 2N. Some of them correspond 
to the intersections of line k y = const with Fermi-surface and have small 
imaginary part of k x ^ a , the others correspond to the point with (Rek x ^ a , k y ) 
far from Fermi-surface. Therefore we obtain from (jSJ) the system of 2iV liner 
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equations for constants C a with E as a parameter. Then the equality of the 
determinant of the system to zero is the condition for existence of bound 
states with energy E: 



( 



U 2 N 



\ 



Det 



U\C 



ik Xl id 



U 2 NG 



V2N 



0. 



(6) 



,ik x<1 {N-X)d 



V2Ne 



ik Xl 2N{N—l)d 



I 



We only consider now the possibility for dispersionless states with E — 
to exist in some region of k y . Then all the solutions of (j3J) have the form 
(u a (k y ) , v a (k y )) T = p a ) with p a = ±1. From each point of inter- 

section k y = const with Fermi-surface we obtain one solution with p a = 
sgn(vf^ x (k x j, k y )A(k x j, k y )) in quasiclassical approximation. And from each 
point far from Fermi-surface we obtain two solutions with close values of k x 
and with opposite values of p a . 

Let n and m be numbers of solutions corresponding to p a = ±1, respec- 
tively. Then we can obtain after some straightforward algebra that in the 
case of n ^ m Eq.© is always true. For n = m Eq.fJUJ) can be reduced to 



The wave vectors k XiCt , corresponding to the same sign of p a , can only 
coincide for a few values of k y , for which different parts or Fermi-surface 
intersect with each other. Therefore we obtain the simple criterion for dis- 
persionless zero-energy bound states to exist: n ^ m. Since the solutions, 
corresponding to the values of k, which are situated far from Fermi-surface, 
always appear in pairs with opposite signs of p a , we can safely take into ac- 
count only solutions with k defined by the intersections of the line k y = const 
with Fermi-surface. 

Let's apply this criterion to the model under consideration. In the qua- 
siclassical approximation we need to obtain all intersections of the line k y = 
const with Fermi-surface in SABZ, and, then, calculate p = sign (vf tX (k) A(k)) 
for all these points. Let us consider all values of k y simultaneously and find 
the positions of the edge of the regions where zero-energy surface states exist. 

Due to the symmetry of the normal metal quasiparticle energy spectrum 
and the superconducting gap to the inversion: £(k) = £(— k), A(k) = A(— k), 



n 






p=-l,i<j 



p=+l,i<j 
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we need to consider only that points at Fermi-surface, where the sign of p 
changes. They are the points of the gap sign changing and the points of v x 
sign changing. 

It is easy to show that points of v x sign changing can't modify the pa- 
rameters n and m. These points are the tangent points of Fermi-surface and 
line k y = const. On the one side (along k y -ax.is) from the point of v x sign 
changing there are two solutions with opposite signs of p a . On the other side 
from this point there are no real solutions, but there are two solutions, which 
have large imaginary part of k x and opposite signs of p a , too. 

Thus let us consider only points of gap sign changing. We should take 
into account univocal correspondence the SABZ and the usual first BZ. 
Then there are only four points of gap sign changing in SABZ, just as in 
BZ: (k®, k®) = (±z/7r, ±i/7r) in crystal axes. Here parameter v takes value 
7T _1 arccos(— fJ,/4t) for the nearest neighbors model and 

v = 7T _1 arccos f ^ | (8) 

\2(t+y/W^¥il)J 

for more general case of the next nearest neighbors model. Parameter v is a 
relative coordinate of BZ point, where the gap sign changing takes place and 
correlated with the filling of the band. Maximal and minimal values of v are 
and 1. For the simplest case of half-filling (/i = 0) we get v — 1/2. But v 
is not strictly the filling of the band. 

Now we need to obtain the coordinates ky of the gap sign changing points 
in SABZ. For (/i/0)-orientation: 

then for /^-coordinates of four gap sign changing points 

k° y = {±h ± l)irud. (10) 

Since k y is a crystal momentum, one can move the /^-coordinates of these 
points into the SABZ. Finally, we obtain the following regions of k y where 
ZES exist: 

I I ^ ( kmin i k max ) , ( 1 1 ) 

where 

k mm = mm{\F[(h - l)nud}\, \F[(h + l)nud]\}. (12) 



5 



k max = max{\F[{h - l)nud]\, \F[{h + l)itvd\\}. (13) 

Here F[...] is a function, which shifts argument to the permissible for SABZ 
value: F[k y ] = ({(1/2) + (k y /2nd)} - (1/2)) 2nd, {...}-is a fractional part of 
the argument. From Eq.fjllj) we can see that the region of k y , where zero- 
energy bound states take place, always exists, except for the case of 

k-min k max - (14) 

It is easy to obtain from Eqs. (|llllT3*j) the regions of ZES existence for 
any cases under consideration. In FigJT] results for (210), (310), and (320) 
surfaces are shown. It is important to note, that regions with ZES and 
regions without them are separated by the lines of zero gap (for this values 
of k y superconducting gap vanishes for one of the quasiparticle trajectories, 
forming the state). 




Figure 1: Grey - regions of ZES existence in (k y ,u) plane. Black lines 
are the lines of zero gap. Dashed line corresponds to half-filling fi = 0. 
Parameter v = 7r _1 arccos(— /i/4t) for nearest neighbors model and v = 
7r _1 arccos(— /i/2(t + \Jt 2 — t'/i)) for next nearest neighbors model. 



For half-filled (/?iO)-surface the result can be formulated in the general 
form: zero-energy Andreev surface states appear only for the case of odd h 
and /. Moreover, it is seen from that for any surface orientation one can 
find set of values of u, for which there are no zero energy states. From ()14|) 
we have h + / values of v (and the same number of band fillings): 

u = i/l, j/h, i = l,...,h-l, j = !,...,/-!. (15) 
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Conclusions. The conditions for zero-energy Andreev surface bound state 
to exist are studied for the lattice model of d-wave superconductor. Arbi- 
trary surface orientation is considered for nearest neighbors as well as for next 
nearest neighbors models. The result is very sensitive to the surface orienta- 
tion, and doesn't change continuously under surface-to-crystal angle rotation. 
In particular, for half-filled (ZiZO)-surface zero-energy Andreev surface states 
only appear under the condition that h and / are odd simultaneously. 
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